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Abstract 

The most general structure of an elastic partial wave amplitude when the unphysical 
cuts are neglected is deduced in terms of the N/D method. This result is then matched to 
lowest order, 0(p 2 ), Chiral Perturbation Theory(xPT) and to the exchange (consistent 
with chiral symmetry) of resonances in the s-channcl. The extension of the method 
to coupled channels is also given. Making use of the former formalism, the irir and 
•KV (1=1/2) P-wave scattering amplitudes are described without free parameters when 
taking into account relations coming from the l/N c expansion and unitarity. Next, the 
scalar sector is studied and good agreement with experiment up to ^/s = 1.4 GeV is 
found. It is observed that the ao(980), a and k(900) resonances are meson-meson states 
originating from the unitarization of the 0{p 2 ) %PT amplitudes. On the other hand, 
the /o(980) is a combination of a strong S-wave meson-meson unitarity effect and of 
a preexisting singlet resonance with a mass around 1 GeV. We have also studied the 
size of the contributions of the unphysical cuts to the 7T7t(I=0) and -Ktt(I=1/2) elastic 
S-wave amplitudes from x?T and the exchange of resonances in crossed channels up to 
y/s w 800 MeV. The loops are calculated as in %PT at next to leading order. We find a 
small correction from the unphysical cuts to our calculated partial waves. 
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1 Introduction 



The understanding of the scalar meson-meson strong interaction is still so controversial that 
there is not even a consensus about how many low lying states (with masses <1 GeV) are 
there. The main difficulties which appear in this sector are: First, the possible presence of 
large width resonances, as the /q (400 — 1200) = a in nir scattering or the Kq (900) = k in the 
1=1/2 Kir amplitude, which cannot be easily distinguished from background contributions. 
Second, the existence of some resonances which appear just in the opening of an important 
channel with which they couple strongly, as for example the /o(980) or the ao(980) with 
the KK threshold around 1 GeV. All these aspects make that, for instance, it is not clear 
how many states are present, which is their nature and why simple parameterization of the 
scalar physical amplitudes in terms of standard Breit-Wigner resonances are not adequate, 
as stressed in several works [jl], ||]. 

The conflictive situation for the scalar sector contrasts with the much better understood 
vector channels. In this latter case, one can achieve a sound understanding of the physics 
involved just from first principles Q, namely, chiral symmetry, unitarity and relations com- 
ing from the QCD limit of infinity numbers of colors (Large N c QCD). This is accomplished 
thanks to the leading role of the /?(770) and if* (890) resonances, in accordance with vector 
meson dominance, very well established from particle and nuclear physics phenomenology. 
The issue is whether such a basic understanding for the scalar sector is possible, and, at 
the same time, is able to reproduce the associated phenomenology. Connected with the 
former, it should be also interesting to see if some kind of scalar meson dominance remains, 
in analogy with the above mentioned vector meson dominance. 

There have been many studies of the scalar sector but none in the basic lines we outlined 
before. These studies have led to a variety of models dealing with the scalar meson-meson 
interaction and its associated low energy spectroscopy. The low energy scalar states have 
been ascribed || to conventional qq mesons [||, §], q 2 q 2 states KK molecules [||, 10 1, 

glueballs [11| and/or hybrids [12|. 

One can think in two possible ways in order to avoid the former explicit model depen- 
dence for the scalar sector, apart, of course, of solving QCD in four dimensions for low 
energies which, nowadays, is not affordable. 

One way to proceed is to make use only of general principles that the physical amplitudes 
must fulfill, as unitarity and analyticity. There are a series of works by Pennington, Morgan 
et al. [|l], |l3| which fit nicely the experimental data available for the scalar-isoscalar sector 
and try to obtain also some understanding of the associated spectroscopy. Another work in 
this line is presented in [fT4[| . However, these approaches have also problems as, for example, 
the specific way in which the amplitudes are parametrized and the lack of enough precision 
in the experimental data in order to discard other possible solutions. 

Another alternative is the use of effective field theories which embody and exploit the 
symmetries of the underlying dynamics, in this case QCD. In this sense, Chiral Perturbation 
Theory (xPT) [jD^, 16 1, is the effective field theory of QCD with the lightest three quark 
flavors. This approach has been extensively used in the last years for the meson sector and 
allows to calculate any physical amplitude in a systematic power momentum expansion. 

This latter point of view will be the one adopted here. First, we will derive, making use 
of the N/D method [17|, the most general structure for an arbitrary partial wave amplitude 
when the unphysical cuts are neglected. In this way, our method can be seen as the zero 
order approach to a partial wave when treating the unphysical cuts in a perturbative sense. 
We think that this will be the case, at least, in those partial wave amplitudes which are 
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dominated by unitarity and the presence of resonances in the s-channel with the same 
quantum numbers of the partial wave amplitude. For example, the case of the p and K* 
resonances in the P wave tttt and Kir scattering respectively or the scalar channels with 
isospin 0, 1/2 or 1, where several resonances appear. In fact, for the p(K*) meson- meson 
channels, at least up to y/s = 1.2 GeV, one can describe accurately the associated tttt^Kit) 
phase shifts just in terms of simple Breit-Wigner parameterizations with the coupling of 
the p(K*) with Tnr(Kn) given by the KSFR relation [18] and their masses taken directly by 
experiment. In this way one has a free parameter description for these processes without 
the unphysical cuts, since it only has the physical or right hand one as required by unitarity. 
Thus, one can deduce that for these processes the contribution from the unphysical cuts is 
certainly much smaller than the one coming from the exchange of these resonances in the 
s-channel and from unitarity. Otherwise, a free parameter reproduction of these channels 
with only the right hand cut could not be possible. This type of description for the p and K* 
meson-meson channels is given below in section 3 and it was also given in ref. 0J for the case 
of the p. For the scalar channels with 1=0,1 and 1/2 there have been a number of previous 
studies ^ ^, ||, IS, 2C] which neglect the contribution from the unphysical cuts establishing 



clearly the great importance of unitarity for the scalar sector. In particular, in the work 
of ref. [19] the 1=0,1 S-wave channels were described in terms of just one free parameter 
up to y/s ~ 1.2 GeV, indicating that the contribution from the unphysical cuts should be 
small enough to be reabsorbed in this free parameter. The connection between the work of 
ref. (l^] and the present one is discussed in section 4 before the subsection dedicated to 
the resonance content of our S-wave amplitudes. It is also interesting to indicate that in 
[21, 22] the S-wave scattering was also studied including the unphysical cuts up to 0(p 4 ) 
in xPT and the results obtained were very similar to the ones of the former works, refs. 
[19, pC|] , without any unphysical cuts at all. Apart from these considerations, we approach 
in the last section of the present work the influence of the unphysical cuts from xPT and 
the exchange of resonances in crossed channels taking into account the results of ref. [23]. 
In this reference, the tttt and Ktt elastic amplitudes are calculated up to one loop including 
explicit resonance fields pi. In this way the range of applicability of xPT is extended up 
to y/s ~ 700 — 800 MeV. The loops are calculated as in XPT at 0(p 4 ). We conclude that 
the contributions of the unphysical cuts are small and soft enough to be reabsorbed in our 
free parameters in a convergent way when treating the left hand cut in a perturbative way. 
It is important to indicate that such a small value for the influence of the unphysical cuts, 
as can be see in Table 3, is due to a cancellation between the contributions coming from 
the loops and the exchange of resonances in crossed channels. 

After neglecting the unphysical cuts we then match the general structure we obtain from 
the N/D method with the lowest order xPT Lagrangian, 0(p 2 ) fl6|| , and its extension to 
include heavier meson states with spin < 1 |24| ], beyond the lightest pseudoscalars (tt, K, rj). 
Making use of this final formalism, we will study the scalar sector, being able to reproduce 
the experimental data up to about y/s < 1.4 GeV, with s the Mandelstam variable cor- 
responding to the square of the total momentum of the pair of mesons. In order to say 
something for higher energies, more channels, apart from the ones taken here, should be 
added. This is not considered in this work, although it can be done in a straightforward way, 
albeit cumbersome, in terms of the present formalism. We also study the vector tttt and 
Kn (1=1/ 2) scattering and compare it with the scalar sector to illustrate some important 
differences between both cases. 

The main conclusion of the work is that one can obtain a rather accurate description 
of the scalar sector compared to experiment, in a way consistent with xPT and Large N c 



3 



QCD |pE|| , if the tree level structures coming from Large N c QCD for the meson-meson 
scattering are introduced in a way consistent with xPT and then are properly unitarized 
in the way we show here. Contrary to what happens for the vector channels, where the 
tree level contributions determine the states which appear in the scattering, we will see 
that for the scalar sector the unitarization of the 0(p 2 ) xPT amplitude is strong enough 
to produce meson-meson states, as for example, the cr(500), ao(980), k(900) and a strong 
contribution to the /o(980). All these states, except the /o(980), disappear for Large N c 
QCD because they originate from effects which are subleading in \/N c counting rules (loops 
in the s-channel). We will see below that the origin for such a different behavior between 
the p and the a will be just a numerical factor 1/6 between the P and S-wave 0(p 2 ) xPT 
amplitude. Note that in a n-loop calculation this gives rise to a relative suppression factor 
of l/6 n+1 of the P-wave loops respect to the S-wave ones. 



2 Formalism 

Let us consider in the first place the elastic case, corresponding to the scattering of two 
particles of masses mi and ra 2 respectively. We will also allow for several coupled channels 
at the end of the section. 

Let T£(s) be a partial wave amplitude with isospin I and angular momentum L. Since 
we are dealing with tt, K and rj as asymptotic particles, which have zero spin, L will also be 
the total spin of the partial wave. The projection in a definite angular momentum is given 
by: 

T[(s) = j\co S dT I ( S ,co S d)P L (cose) (1) 

where (V2) a is a symmetry factor to take care of the presence of identical particles states 
as rjrj or tttt, this last in the isospin limit. The index a can be 0,1 or 2 depending of the 
number of times these identical particle states appear in the corresponding partial wave 
amplitude. For instance, a = 2 for tttt — > tttt, a = 1 for r/r/ — > KK, a = for Ktt — > Ktt 
and so on. Pl(cos0) is the Legendre polynomial of L th degree. 

A T[(s) partial wave amplitude has two kinds of cuts. The right hand cut required by 
unitarity and the unphysical cuts from crossing symmetry. In our chosen normalization, the 
right hand cut leads to the equation (in the following discussions we omit the superindex 
/, although, it should be kept in mind that we always refer to a definite isospin): 



ImT^ 1 = -p(s) (2) 

for s > sthreshoid = $th- In the case of two particle scattering, the one we are concerned 
about, s t h = (wi + rrt'i) 1 and p(s) is given by: 

Pis) = (3) 

with 



_ V(s - (mi + m 2 ) 2 )(s - (mi - m 2 ) T ) = A 1/2 (s, mf, m|) 
P ~ 2^ 2^5 U 

the center mass (CM) three momentum of the two meson system. 
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The unphysical cuts comprise two types of cuts in the complex s-plane. For processes 
of the type a + a — > a + a with mi = m-2 = m a , there is only a left hand cut for s < SLeft- 
However for those ones of the type a + b — > a + b with mi = m a and m,2 = m&, apart from 
a left hand cut there is also a circular cut in the complex s-plane for \s\ = m| — mf, where 
we have taken m-2 > mi- In the rest of this section, for simplicity in the formalism, we will 
just refer to the left hand cut as if it was the full unphysical cuts. This will be enough for 
our purposes in this section. In any case, if we worked in the complex p 2 -plane all the cuts 
will be linear cuts and then only the left hand cut will appear in this variable. 

The left hand cut, for s < SLeft, reads: 

T L (s + ie)-T L (s-ie) = 2ilm T L (s) (5) 

The standard way of solving eqs. (||) and (||) is the N/D method flT^ ]. In this method 
a T^(s) partial wave is expressed as a quotient of two functions, 

Ti(s) = BUT) (6) 

with the denominator function Dl(s), bearing the right hand cut and the numerator function 
Nx,(s), the unphysical cuts. 

In order to take explicitly into account the behavior of a partial wave amplitude near 
threshold, which vanishes like p 2L = v L , we consider the new quantity, T' L , given by: 

W = ^r 1 (7) 

which also satisfy relations of the type of eqs. (||) and (|5|). So that we can write: 

r ( S ) - n ' l{s) (8) 
Tl[s) ~ nfisj (8) 

From eqs. (||), (|5|) and ([/]), N^(s) and D' L (s) will obey the following equations: 
Im D' L = Im T' L l N' L = -p{ S )^' L v\ s > s th 

Im B' L = 0, s < s th 1 } 

Im N' L = Im T' L B' L , s < s Left 

ImN' L = 0, s> s Left 1 j 

Since and D' L can be simultaneously multiplied by any arbitrary real analytic function 
without changing its ratio, T' L , nor eqs. (^) and (10), we will consider in the following that 



is free of poles and thus, the poles of a partial wave amplitude will correspond to the 
zeros of D' L . 

Using dispersion relations for D^(s) and N^(s), we write from eqs. (|) and ©: 



n— 1 



where n is the number of subtractions needed such that 

hm ^4 = (12) 
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since, from eq. (jH) 



lim 



v L p(s) 



1 



4 L + 2 ir 



On the other hand, from eqs. and (0), consistently with eq. (12) 

\n-L t m / t\ t~v/ / A n.-L-l 



N' L ( S ) 



(s - so) r 



7T 



Im T L (s')D^') 



v(s') L (s' - s ) n - L (s' - s) 



+ E 



(13) 



(14) 



m=0 



Eqs. (|ll]) and ( |14] ) constitute a system of integral equations the input of which is given 
by Im Tl(s) along the left hand cut. 



However, eqs. (11) and (14) are not the most general solution to eqs. (B) and ( |10| ) 
because of the possible presence of zeros of Tl which do not originate when solving those 
equations. These zeros have to be included explicitly and we choose to include them through 
poles in D' L (CDD poles after ref. |26|| ). Following this last reference, let us write along the 
real axis 



Im B' L (s) 



d\(s) 
ds 



(15) 



Then by eq. © 



— = -p(s)U L W L , S > 8th 



dX ~ () 

ds 



(16) 



s < S t h 



Let Si be the points along the real axis where T^(sj) = 0. Between two consecutive 
points, Si and Sj+i, we will have from eq. fll(i| ) 



AW 



u(s') L p(s')^ L (s')ds' + X( Si ) 



(17) 



with A(sj) unknown because the inverse of T' L (si) is not defined. Thus, we may write 



X(s) 



I v(s') L p(s')N' L (s')ds' + ^A( Si )0(s 



(18) 



with 9(s) the usual Heaviside function. 
From eqs. (P) and (pi), it follows that 



(s-so) n f°° Im D / L (s / )dg / 
„ ( s '- s )( 5 '- So ) 



7T 



71-1 



n + E 



— m 

(Xrn.S 



m=0 



n-1 

£ 

m=0 



— m 

u, m S 



7T 



th (s'-s)(s'- So ) n 



(s _ SQ )n Z" 00 Z/(s')^/3(s')N^(s') , 



7T 



7T 

n-1 



(s'-s)(s'-s )' 



(s'-s)(s'-s ) r 



^'+E 



a m s 



m=0 



+ E 



Afa) (g - s ) n 
r' 7r(si - s ) n Si - s 



(19) 
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Eq. (19) can also be obtained from eq. (g) and use of the Cauchy theorem for complex 
integration and allowing for the presence of poles of (zeros of T^) inside and along the 
integration contour, which is given by a circle in the infinity deformed to engulf the real axis 
along the right hand cut, sth < s' < oo. In this way one can also consider the possibility of 
there being higher order zeros and that some of the Sj could have non-vanishing imaginary 
part (because of the Schwartz theorem, s* will be another zero of T' L (s)). However, as we 
will see below for L < 1, when considering chiral symmetry in the Large N c QCD limit, 
the zeros will appear on the real axis and also as simple zeros. In general, using T' L instead 
of Tl, we avoid working with L th order poles of at threshold in the dispersion relation 
given by eq.(|i~9|). 

The last term in the right hand side of eq. (|i~9| ) can also be written in a more convenient 
way avoiding the presence of the subtraction point so- To see this, note that 

(S ~ S ) n -tS-Si+Si-So 

— — = (s-so) — 

S Si S S~i 

, \n-l/i , s i~ s 0\ I \n-\ i 1 \( s ~ s o) n 

= {s-sq) (H — — ) = (s-s ) +{si-s ) — - — 

S Si S Si 

n-l 



i=0 



s-s r- l - i (s i -sof + ^ S ° r (20) 

S - Si 



The terms 



can be reabsorbed in 



As a result we can write 



n-l 
i=0 



n-l 

E' 

m=0 



(»-»„)" r "W) l pww l (s') , ^i- 7i ,,,, 
d^w- — — y w(s ,_ s)(s ,_ So) „ + E«".« +Ej^- PD 

with a m (n — 1 > m > 0) and 7,, Sj (i > 0), arbitrary parameters. However, if some of the 
Si is complex there will be another Sj such that Sj=s* and jj = 7*, as we explained above. 
Each term of the last sum in eq. (^) is referred as a CDD pole after [^] . 

Eqs. ( |21]) and ( p^ ) stand for the general integral equations for D' L and N' L , respectively. 
Next we make the approximation of neglecting the left hand cut, that is, we set Im T^(s) = 
in eq. (|T4]). Thus one has: 

n-L-l 

Ntw = E s - sm ( 22 ) 

m=0 

As a result, N' L (s) is just a polynomial of degree < n — L — 1 Q . So we can write, 

n-L-l 

N^( 5 )=C J] ( 23 ) 
i=i 

1 One can always make that n > L + 1 just by multiplying and D' L by s fc with A; large enough. 
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In eq. ( p3|) it is understood that if n — L — 1 is zero is just a constant. Thus, the 
only effect of will be, apart of the normalization constant C, the inclusion, at most, of 
n — L — 1 zeros to T' L (s). But we can always divide N' L and by eq. (^). The net result 
is that, when the left hand cut is neglected, it is always possible to take N^(s) = 1 and all 
the zeros of T^(s) will be CDD poles of the denominator function. In this way, 
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The number of free parameters present in eq. (|24| ) is L + l + 2^, where g is the number of 
CDD poles, Ml, minus the number of complex conjugate pairs of Sj. These free parameters 
have a clear physical interpretation. Consider first the term 2g which comes from the 



presence of CDD poles in D' L (s), eq. (24). In [27] the presence of CDD poles was linked 
to the possibility of there being elementary particles with the same quantum numbers as 
those of the partial wave amplitude, that is, particles which are not originated from a given 
'potential' or exchange forces between the scattering states. One can think that given a 
~D' L (s) we can add a CDD pole and adjust its two parameters in order to get a zero of the real 
part of the new D^(s) with the right position and residue, having a resonance/bound state 
with the desired mass and coupling. In this way, the arbitrary parameters that come with 
a CDD pole can be related with the coupling constant and mass of the resulting particle. 
This is one possible interpretation of the presence of CDD poles. However, as we are going 
to see below, these poles can also enter just to ensure the presence of zeros required by 
the underlying theory, in this case QCD, as the Adler zeros for the S-wave meson- meson 
interaction. The derivative of the partial wave amplitude at the zero will fix the other CDD 
parameter, 7^. With respect to the contribution L + 1 to the number of free parameters 
coming from the angular momentum L, it appears just because we have explicitly established 
the behavior of a partial wave amplitude close to threshold, vanishing as v L '. This is required 
by the centrifugal barrier effect, well known from Quantum Mechanics. 

It should be stressed that eq. ( f24| ) is the most general structure that an elastic partial 
wave amplitude, with arbitrary L, has when the left hand cut is neglected. The free param- 
eters that appear there are fitted to the experiment or calculated from the basic underlying 
theory. In our case the basic dynamics is expected to be QCD, but eq. (p4j) could also 
be applied to other scenarios beyond QCD as the Electroweak Symmetry Breaking Sector 
(which also has the symmetries [|28| used to derived eq. (p4[), as far as it is known). 

Let us come back to QCD and split the subtraction constants a m of eq. in two 

pieces 

a m = ai + ai L {s ) (25) 

The term will go as N c , because in the N c — > 00 limit, the meson-meson amplitudes go 
as iV" 1 [p5fl . Since the integral in eq. (|24] ) is 0(1) in this counting, the subleading term 
a m L ( s o) is of the same order and depends on the subtraction point so- This implies that eq. 
(|24|), when N c — > 00, will become 
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B' L (s) = (,) = 53 a£s m + ^ -i5L_ (26) 

S — Si 

m=0 i 

where Rf is the leading part of Ri and M£° counts the number of leading CDD poles. 



Clearly eq. (26) represents tree level structures, contact and pole terms, which have 
nothing to do with any kind of potential scattering, which in Large N c QCD is suppressed. 

In order to determine eq. ( p6| ) we will make use of xPT [pj and of the paper |24| j. In 
this latter work it is shown the way to include resonances with spin < 1 consistent with 
chiral symmetry at lowest order in the chiral power counting. It is also seen that, when 
integrating out the resonance fields, the contributions of the exchange of these resonances 
essentially saturate the next to leading xPT Lagrangian. We will make use of this result 



in order to state that in the inverse of eq. (26) the contact terms come just from the lowest 



order xPT Lagrangian and the pole terms from the exchange of resonances in the s-channel 
in the way given by |Q| (consistently with our approximation of neglecting the left hand cut 
the exchange of resonances in crossed channels is not considered) . In the latter statement it 



is assumed that the result of [24] at 0(p A ) is also applicable to higher orders. That is, local 
terms appearing in xPP anci from eq. ( |26| ) of order higher than C(p 4 ) are also saturated 
from the exchange of resonances for N c >> 1, where loops are suppressed. 

Let us prove that eq. (|2q) can accommodate the tree level amplitudes coming from 



lowest order xPP [ 16 1 and the Lagrangian given in |24| for the coupling of resonances (with 
spin < 1) with the lightest pseudoscalars (n, K and rf). 

Following p4|] , one can write the exchange of a resonance (with spin < 1) divided by v L 

as: 



2\2 



(ciS + c-m 2 ) 
di.s 



M 2 



L = 



L = 1 (27) 



where Mj is the mass of the i th resonance, m 2 is some combination of squared masses of 
the lightest pseudoscalars and q, d i and di are arbitrary constants with di > 0. 

The lowest order xPT partial wave amplitude (L < 1), divided again by u L , can be 
written schematically as: 

as + a'rh 2 ; L = 

b ; L = 1 (28) 

with rh? another combination of squared masses of the lightest pseudoscalars. 

Thus, in Large iV c QCD, the partial wave amplitudes will have the following structure, 
after omitting the exchange of resonances in crossed channels and contact terms of order 
higher than 0(p 2 ) 



T o( s ) = T o (s)=as + am + ^ — M 2 _ s 



i=i 

%{s)^r{s) = b + i^ 1 ^- s (29) 

i=i 1 
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In the former equation it is understood that if Rl = the sum does not appear. 
Let us now study the inverse of eq. ( |29|) in order to connect with D^ . With the 
number of zeros of T^°, we can write: 

i A nt\(s-M?) 
^ris)- AL uius-s r ) (30) 

In the former equation, if Rl = or ^ = 0, the corresponding product must be 
substituted by 1. Al is just a constant. Note that from eqs. (]29|) £o < Ro + 1 and £i < 
from simple counting of the degree of the polynomials that appear in the numerators after 
writing eqs. (^) as rational functions. 

On the other hand, note from eqs. ( |27| ) that, for s > 0, the amplitude from the exchange 
of resonances, both for L = and 1, is positive below the mass of the resonance and negative 
above it. This implies that in the interval, Mf, 1 > s > Mf, by continuity, there will be a 
zero of T^°. In this way, 

Zl > Rl - 1 (31) 

For L = 0, apart from the zeros between resonances, one has also the requirement from 
chiral symmetry of the presence of the Adler zero, along the real axis and below threshold. 
Thus, £o > Ro and we can write 

Rl-L + 1>£ l >R l -L (32) 

Note that since the Rl — L zeros of are real, the possible single additional zero from 
the upper limit of eq. (32) must be also real, since a complex one would imply its conjugate 



too. This is so because all the coefficients in eq. (2£) are real. 

Let us do the counting of zeros and resonances from eq. (|26|). The number of zeros of 
T£° is equal to the number of CDD poles of D£°, Mf. Hence 

£l = M£° (33) 

On the other hand the number of poles of T'£° is equal to the number of zeros of T)'£° 
which has the following upper limit 

Rl<MF + L = Z l + L (34) 



^l>Rl-L (35) 



Thus, 



which gives the same lower limit as in eq. (|32[). Let us recall that the zeros of eq. ( p0|) 
are real, as discussed above. This, together with the requirement of T^°(s) being a real 
function on the real axis, forces all parameters in eq. (^) to be real. Then the number of 



free parameters in eq. ( 26 ) are 2M£° + L + 1 = 2£l + L + l. When fixing the parameters 
in eq. (^) to the s r parameters in eq. ([}(]) this reduces in £x the number of free parameters 
in eq. (|^). By fixing the arbitrary constant Al of eq. ( |30| ) this leaves us with £l + L free 
parameters in eq. (p6|). Imposing now the position of the Rl resonances of eq. ( p0| ) to be 
at Mf we have Rl additional constraints. Hence we should have 

£ L + L - R L > (36) 
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which actually holds, as seen above in eq. (|3~5l). As a consequence eqs. ( ^S|) can always be 



cast in the form of eq. (26) 



Let us define the function #l(s) by 

\L+1 roo r,(J\L, 



g L (s)u -^a m (s )s ^ ds _ _ Sq)L+1 (37) 



m=0 Js th 

and using the notation 



L M L 

Tf(s) = aL ^ m + E -^r 1 (38) 

z — ' z — ' S — Si 

m=0 i=l 

results that from eq. (|3]) one has 

T L ( S ) = [l/T2°( S )+ 5i ( S )]- 1 (39) 

The physical meaning of Eq. ([^) is clear. The T^ 3 amplitudes correspond to the 
tree levels structures present before unitarization. The unitarization is then accomplished 
through the function It is interesting at this point to connect with the most popular 

K-matrix formalism to obtain unitarized amplitudes. In this case one writes: 

T L (s) = [K L {s)- 1 - ip(s)}- 1 (40) 

TOO —1 / 



we see that the former equation is analogous to eq. (j3S| ) with T^ 1 ~ (s)+Re (?l(s) = Kl{s) . 

From the former comments it should be obvious the generalization of eq. ( p9| ) to cou- 
pled channels. In this case, T^(s) is a matrix determined by the tree level partial wave 
amplitudes given by the lowest order %PT Lagrangian [|0| and the exchange of resonances 
||. For instance, [T£°(s)]n = T$ +T$, [T£°(s)]i 2 = T$ and so on. Where is 

the matrix of the 0(p 2 ), xPT partial wave amplitudes [10]] and T R the corresponding one 



from the exchange of resonances [24] in the s-channel. Once we have T^s) its inverse is 
the one which enters in eq. (|39|). Because N^(s) is proportional to the identity, ^(s) will 
be a diagonal matrix, accounting for the right hand cut, as in the elastic case. In this way, 
unitarity, which in coupled channels reads (above the thresholds of the channels % and j) 

[Im T L l ]ij = -pn{s)5ij = Im g L (s)u5ij (41) 

is fulfilled. The matrix element gL{s)u obeys eq. (j37|) with the right masses corresponding 
to the channel i and its own subtraction constants af L (so). 

In the present work the coupling constants and resonance masses contained in T^(s) 
are fitted to the experiment. The same happens with the af L although, as we will discuss 
below, they are related by SU{3) considerations. 

In Appendix A the already stated coupled channel version of eq. (39) is deduced directly 
from the N/D method in coupled channels ]p^| . 



3 The ideal case: elastic vector channels. 

In this section we are going to study the tttt and Kir scattering with I=L=1 and 1=1/2, L=l, 
respectively. It will be shown, as mentioned in the introduction, that these reactions can 
be understood just by invoking: 
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1) Chiral Symmetry 

2) Large A^ C QCD 

3) Unitarity 

A priori, one could think that relations coming from the N c — > oo limit should work 
rather accurately at the phenomenological levels in these channels due to the predominant 
role of the p and K* poles over subleading effects in 1/N C , as unitarity loops. 

In eq. ( f24| ) the L=l zero at threshold is included in the Ni function, eq. @, which can 
be taken as 



Ni 
Di 



i/Ni 



(42) 



with and Dj defined in eq. ® and given in eq. (24). However, in this chapter we will 
also include the P-wave threshold factor as a CDD pole in the Di function. In this way, we 
will treat the S-wave Adler zero and the P-wave threshold zero in the same footing, which 
will make the comparision among both partial waves more straightforward. The final result 
can be derived in the same way as eq. (24) but working directly with Tj instead of T' x . 
Thus: 



Ni(s) 
DiOO 



(43) 



E 



+ a 



so 



Si 



7T 



p(s') 



ds'- 

(s'-s)(s'-s ) 



Of course, the partial wave amplitude Ti will be the same than before since we have 
just divided at the same time Ni and Di by v = p 2 . As a simple and explicit example that 



eq. fl24|) together with eq. ([42|) are equivalent to eq. (|43|), let us consider the scattering 
of two pions. After multiplying Ni and Di in eq. fl42] ) by 4, we divide both functions by 
4 v = s — sa, with sa = 4 m 2 . Keeping in mind eq. (^4|), one still has a sum over the former 
CDD poles together with another one at s = sa- The subtraction polynomial of order one 
in s transforms into a CDD pole at s = sa plus a constant. Finally, the dispersive integral 



gives rise to that of the former equation (43) together with a constant plus a CDD pole at 
s = sa — to see this just add and subtract sq in 4i/(s') in the integral of eq. (|24"|). Hence, 



the structure given in eq. (42) is obtained 



q 

P / m i P 

P-c, ^ 

Figure 1: Loop giving rise to the go function. 

The integral in eq. (^) will be evaluated making use of dimensional regularization. It 
can be identified up to a constant to the loop represented in Fig. 1. This identification is 
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consequence of the fact that both the integral in eq. (^) and the loop given in Fig. 1 have 
the same cut and the same imaginary part along this cut, as it can be easily checked. 
Following eq. ([37]) we define 



90 0) 



a SL (so) 



so 



7T 



p(s') 



ds'- 

!h (s> ~ s)(s> - so) 



(47T) 



•log 



~a SL 



(ji) + log -f 



m 



to 2 + s m\ X 1 ^ 2 (s,ml,ml) 



2s 



log -4 



2s 



m\ + m\ 



s + A 1 / 2 (s, ml, ml) 



m 2 + to 2 . 



X 1 / 2 (s, mf , m§) 



(44) 



for s > s^/j. For s < or s complex one has the analytic continuation of eq. (p4|). The 
function A 1 / 2 (s, mf , to 2 ,) was already introduced in eq. @. The regularization scale 
appearing in the last formula of eq. (K3|) , plays a similar role than the arbitrary subtraction 
point so i n the first formula of eq. (|44l). This similarity is consequence of the fact that both 
\i and so can have any arbitrary value but the resulting function g§ (s) is independent of this 
particular value because of the change in the subtraction constant, a SL (/j.) for dimensional 
regularization or a SL (so) for the dispersion integral. The a SL (n) 'constant' will change 
under a variation of the scale \i to other one \j! as 



~SL 



~a SL 



{H) + log ^ 
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(45) 



in order to have go( s ) invariant under changes of the regularization scale. We will take 



/i = Mp = 770 MeV |^0[. The function go(s) is also symmetric under the exchange mi m2 
and for the equal mass limit it reduces to 

1 



9o(s) 



(47T) 



a 5L (/i)+log^-+fj(s)log 



I 1 



a(s) + 1 
a(s) - 1 



with 



a(s) 



1 



4m 2 



(46) 



(47) 



After this preamble, let us consider in first place the tttt scattering with I=L=1. As can 
be seen from Fig. 2 this process is dominated by the p exchange. From eq. (E3) one has 



7i 



4to 2 



i=2 



Si 



-1 



(48) 



The first term in the R.H.S of the last equation fixes the zero at threshold for a P-wave 
amplitude. 



The tree level part of T^ n (s) from ref. [24] and lowest order X-PT [16], in the way 
explained in the last section, is given by: 



+ q 2 



2 p. 



3 / 2 M 2 - s 



(49) 



with p 2 ,,. the three-momentum squared of the pions in the CM, / = 87.3 MeV the pion 
decay constant in the chiral limit . The deviation of g% with respect to unity measures 
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the variation of the value of the p coupling to two pions with respect to the KSFR relation 
]l8[] , g 2 = 1. In |3l| this KSFR relation is justified making use of Large N c QCD (neglecting 
loop contributions) and an unsubtracted dispersion relation for the pion electromagnetic 
form factor (a QCD inspired high-energy behavior). 

Comparing eqs. ( |48| ) and (^), one need only one additional CDD pole apart from the 
one at threshold and we obtain 



72 



«2 







6/ 2 (M p 2 



4m: 



{Mj-Aml(l-gD) 
6/ 2 9lMl 



1 - g 2 M 2 - (1 - g*)4m% 

M 2 
p 



(50) 



Thus, we can write our final formula for the isovector tttt scattering in the following way 

-l 



7i 



4 m 2 



+ 



'v 7r7r 
72 

3 - S 2 



+ 9o 



(51) 



in terms of the parameters g 2 and a (fi). Since g 2 is expected to be close to unity as 
discussed above, it is useful to consider the limit when g 2 — > 1, in which case s 2 —>■ oo such 
that 



72 



S2 



6^ 

M 2 



(52) 



in which case the second CDD pole in eq. (^), at s 2 , moves to infinity and the CDD pole 
contribution gives rise to a constant term, a'^. In this limit we can write 



to 



4m 2 



(53) 



Our calculated phase shifts for the vector tttt scattering are represented in the dashed 
line of Fig. 2, when g 2 is taken equal to one and a SL = is set to zero. The agreement with 
the experimental data is rather good. If we had just taken the imaginary part of gQ W (s) 
in the former equation we would have obtained basically the same curve. This last case 
corresponds just to the Breit Wigner amplitude 



2 v 2 



M 2 
p 



3/2 m 2 -s- \M p Y p {s) 
I» = 



(54) 



12tt/ 2 y/Ts 

For the 1=1/2, L=l Ktt scattering, the tree level amplitude Tf^ 00 , is just given by 

777*^ TTi^ ^ 

multiplying eq. (f49|) by 3/4 and substituting p 2 w by p 2 Klv = - — ^ — — and M p by 



4s 



M^-=896 MeV pOfl, the mass of the neutral K* (890). Now m,K is the mass of the kaon, 
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Figure 2: Isovector irir elastic phase shifts from threshold up to < 1.2 GeV. The dashed 
line corresponds to take g 2 = 1 and a 5L = 0. The continuum line corresponds to the 
simultaneous fit to the p and K* channels, given by eq. (|||). Data: circles [p2;1 , squares 



495.7 MeV [3C]. On the other hand, since 7^ m v instead of having just a single pole 
in the denominator function, as in eq. (|48|), for the zero at threshold, one has two simple 
poles 



7 pr 



7 K« 



(55) 

4 PKn 

with s + = (mj- + m^) 2 and s_ = (mj{ — m^) 2 . That is, we will have two CDD poles but 
both entering with just one parameter because the behavior at threshold is proportional to 
In our notation 

Tp(s) 4 



lim 



7p 



For the simple and realistic case, g 2 = 1, we have, analogously to the case of the p, 



+ a^ 7r + ^( S )]- 1 



(56) 



with 



K 71 



Ti 



(57) 
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In the general case when g% ^ 1, one proceeds in the same way as for the p introducing 
an extra CDD, but we shall omit the details here and the evaluations are done directly 
using the final formula 

T ** =[ l/TK*°°+gK*(s) ] 1 (58) 

with T^ noc evaluated as mentioned above. 

In Fig. 3, the calculated phase shifts for the 1=1/2, P-wave Ktt scattering are shown in 
the dashed curve for g% = 1 and 5 = 0. The same remarks, as done before for the p when 
commenting Fig. 2, are also valid for the K* . It is worth stressing that the dashed lines of 
Fig. 2 and 3 have no free parameters at all, and depend only on / and the masses of the 
resonances K* and p. 

200 n 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 n 




700 800 900 1000 1100 



E em [MeV] 



Figure 3: P-wave elastic Ktt phase shifts with 1=1/2, from threshold up to y/s < 1.2 
GeV.The dashed line corresponds to take g% = 1 and a SL = 0. The continuum line corre- 



sponds to the simultaneous fit to the p and K* channels, given by eq. (|59|). Data: triangles 



circles [35] 



The subleading constant a SL presents in go(s), eq. (44), should be the same for the tttt 
and Ktt states because the dependence of the loop represented in Fig. 1 on the masses of 



the intermediate particles is given by eq. (f44|). This point can be used in the opposite sense. 
That is, if it is not possible to obtain a reasonable good fit after setting a SL to be the same 
in both channels, some kind of SU(3) breaking is missing. 

we do a simultaneous fit to the experimental P-wave tttt (1=1) 



From eqs. (51) and ( r 

and Ktt (1=1/2) phase shifts with g% and a as free parameters using the minimization 
program MINUIT. In order to make the data from different experiments consistent between 
each other, a systematic relative error of a 5% is given to each data point if its own error is 
smaller than this bound. The result of the fit is: 
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g 2 v = 0.87883 ± 0.0157 
a SL = 0.34136 ± 0.0422 (59) 

the errors are just statistical and are obtained by increasing in one unit the x 2 P er degree 
of freedom, x\ f ^ ne Xdof obtained is: 

X 2 d . . f = 0.7406 (60) 

with 81 experimental points. 

The continuum line corresponds to the fit given by eq. (|59[). We see that the agreement 
with data is very good. Note also that g% is very close to unity. To consider the discrepancy 
with the KSFR result, which refers to the value of the coupling constant, it is better to use 
g v which results to be, from eq. (|59[), 0.94. That is, only a 6% of deviation respect to unity. 

It is interesting and enlightening to see the value that the regularization scale p should 
have in order to generate the p pole at 770 MeV when removing 5^ from eq. (|5^) and 
setting a SL = 0. In this way we are taking the regularization scale as a cut off. Eq. ((53|), 
neglecting 4m^ with respect to M 2 , transforms to: 



p ' 



Ami 



(61) 



The resulting p will be around by 0.7 TeV, a value completely senseless 0. Its natural 
value is p ~ 1 GeV, where typically resonances appear. A similar conclusion about these 



unrealistic high values of the cut-off was also obtained in |36|]. This value of p makes it 
manifest that, as we have already seen in this section, the origin of the vectors K* and p 
is attached to tree level structures, preexisting before unitarization. This case is different 
than what will be found in the scalar case where unitarity plays a very important role. 

For the scalar sector, which we will study in detail in the next section, we have only to 
make the following change in eq. (pl|), 

6/ 2 _^ / 2 (62) 
s — Am\ s — 

that is, basically a factor 6 of difference for s around M^. This makes that the 'cut-off' 
needed in the irir S-wave to get a resonance of the same mass than the p is just 1.8 GeV. 
The change has been drastic due to the logarithmic dependence of the regularization scale 
in eq. (|6"l"|). We will see below the consequences that follow from this fact for the scalar 
sector. 

From eqs. (48) and (49), one can see the limitations that some unitarization methods 



of the xPT amplitudes can have when handling a general situation. In particular one can 
think of the Inverse Amplitude Method (IAM) |37|. This method has proved to be a very 
powerful tool to extend the range of applicability of the perturbative series of xPT up to 
energies around 1-1.2 GeV, giving rise to the physical resonances that appear in that case 
as the p, and K* in the vector channels as well as the <r, /o(980), k and ao(980) in the scalar 
ones [^, j3S} |. This approach is based in an expansion of the inverse of the K matrix [20], see 

2 Witout neglecting 4mJ compared with M% one has to multiply the quotient s ®[ m i in eq. (^) by 

— 4n 



M —Am I — ~T 

0.87, as one can check from eq. The resulting fi is 0.3 TeV 
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eq. ([!(]) and the sentence below it. Consider that one can neglect the loop contributions in 
the K matrix, as it occurs for the tttt and 1=1/2 Kir P-wave amplitudes and it is also the 
case in large N c QCD. Then, from eq. (p9|), if gl 7^ 1 a zero of the tree level amplitude will 
appear, which, in turn, implies a pole in its inverse. As a consequence, the expansion of the 
inverse amplitude will break at (note that from eq. ( M ) a zero i n the K matrix implies a 
zero in the full amplitude Tx): 

Ml 

s 2 = ^ (63) 

where the zero of the amplitude occurs. 

Thus, for our final value g% ~ 0.9, ^/s2 ~ 2.4 GeV, further away of the physical region 
typically studied by this method. That explains the success of the IAM when applied to 
the vector tttt and 1=1/2 Kir channels in the region dominated by the p and K* resonances, 
respectively. However, in a general scenario with g% very different from 1 we cannot guar- 
antee the applicability of the IAM. This point should be considered when making use of 
the IAM in other situations beyond QCD, as the ESBS pGj] , where the underlying theory 
is not known and, furthermore, there are no experimental data to rely upon. 



4 The scalar sector. 

In this section we want to study the S-wave 1=0,1 and 1/2 amplitudes. For the partial 
wave amplitudes with L=0 and 1=0 and 1, coupled channels are fundamental in order 
to get an appropriate description of the physics involved up to y^s < 1.4 GeV. This is 
an important difference with respect the former vector channels, essentially elastic in the 
considered energy region. Up to ^fs = 1.4 GeV the most important channels are: 

1 = 7rvr(l), KK(2), 7/7/(3) 

I = 1 7r?7(l), KK{2) (64) 
I = 1/2 Ktt(1), Kr)(2) 

where the number between brackets indicates the index associated to the corresponding 
channel when using a matrix notation as the one introduced in section 2. 

For the 1=0 S-wave, the 4-7T state becomes increasingly important at energies above 
1.2 — 1.3 GeV, so that, in this channel, we are at the limit of applicability of only two 
mesons states when -y/i is close to 1.4 GeV. In the 1=1/2 channel, the threshold of the 
important Krj' state is also close to 1.4 GeV. Thus, one can not go higher in energies in a 
realistic description of the scalar sector without including the Kr( and 4ir states. 

Two sets of resonances appear in the former L=0 partial wave amplitudes IpOf . A first 
one, with a mass around 1 GeV, contains the 1=0 /o(400 — 1200) and /o(980) and the 
1=1 ao(980). A second set appears with a mass around 1.4 GeV as the 1=0 /o(1370) and 
the /o(1500), the 1=1 ao(1450) or the 1=1/2 Kq(1430). As a consequence, one could be 
tempted to include the exchange of two scalar nonets, with masses around 1 and 1.4 GeV. 
Before discussing whether this is the case, let us write the symmetric T§° matrix of tree 
level amplitudes for the different isospins. This matrix, T§°, is determined, as explained at 
the end of section 2, from the lowest order xPT amplitudes, T^ 2 \ and from the exchange 
of scalar nonets in the s-channel as given by ref. pi], T^. In the following formulas we 
consider the exchange of only one nonet. If more nonets are needed, they have only to be 
added in the same way as the first nonet is introduced 
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1=0 



rpOO 

X 0,U 



rpOO 

L 0,12 



rp OO 

1 0,22 



rpOO 

L 0,13 



r-pOO 

^0,23 



rpOO 

J-0,33 
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2/2 



2Mf 



2 Mf 
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s Mf 
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^3 



ai«3 



V3/3(0)i/?(0) 3 



12/2 



3s - 2m 2 rj 
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- 2/3m 2 



Ml 
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16m 



K 



4//„ 
7m? 



Ml 



18/ 2 



1 (a 3 ) 2 



2Mf 



+ 1 (/?(0) 3 ) 



2 Ma 



(65) 



with Mi and M$ the masses of the singlet and octet in the SU(3) limit, is the mass 
of the rj, 547.45 MeV [Q, f v is the decay constant of the 77, set to the value = 1.3 f n 
according with the xPT prediction [40 1 , f n = 93.3 MeV is the pion decay constant and ctj 
and P(0)i are given by: 



/3(0)i 
/3(0) 2 

"2 
0.3 



mz 



Cd)m K - (5c m - c d )— -] 



4 s 

72^2 + ^ 



y2 2 



72^2 + ^ 



(66) 



the constants q, c m , and c m characterize the coupling of a given scalar nonet to the 



pseudoscalar pairs of pions, kaons and etas as given in [24|. 



1=1 



rpOO 



rpOO 

1 0,12 



rpOO 

x 0,22 



ml (/?(l)i) 2 
3/ 2 + Ml-s 

s , (l(l) 2 ) 2 
4/2 + 'm 2 - S 



V2 



Ml - s 



(67) 
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with the function /3(l)i given by: 



\/2 



P(lh = ^[cd^ + icm-cjmjc] (68) 



1=1/2 



5s 2 - 2s(ml + m\) - 3(m| - m 2 ) 2 3 (P{1/2) 1 ) 2 



-9s 2 + 2sm 2 K + 3sm 2 + 7sm 2 - 9m\ + 9m 2 K {m 2 x + m 2 ) - 9m\m\ /3 0(l/2)i/3(l/2) 2 



T' 

24,/! "V2 M| 

-9s 2 - 9(m 2 ~ ^') 2 + ^(3m\ + ™g) ~ 4sm| Q3(l/2) 2 ) 2 
T °> 22 " 2l7f| + M|- S (69) 



with fjc the kaon decay constant with the value fx = 1-2/^ according to experiment [41]. 
The functions /3(l/2)j are given by 



1 

2 o_2\ „_/_2 , _2^ 



/3(l/2)i = -^-[c d s + (c m -c d ){m\ + m 2 )] 
/if 

/3(l/2) 2 = -l—[c d s + c m {5m 2 K - 3ml) - c d {m 2 K + m 2 )} (70) 

Note that the introduction of a nonet implies six new parameters, two masses and four 
coupling constants, which we fit to the experiment. 



According to eq. (p9|), we also need the function go(s), given by eqs. (44) and (jig), with 
its corresponding a for the S-wave channels. By SU(3) arguments, the a SL constant can be 
different for vector and scalar channels. The reason is that a two meson state has different 
SU(3) wave functions in S and P wave, because under the exchange of both mesons the 
spatial P-wave is antisymmetric while the S-wave is symmetric and the total wave function 
must be symmetric. That is, the two mesons are in different SU(3) representations. 

We have included /„ in the S-wave isoscalar 0(p 2 ) X-PT amplitude for KK — ► rjrj and 
fx in the S-wave, 1=1/2 tree level amplitudes, to obtain, after the fit, that the a constant 
was the same for all the scalar channels. These changes come from the SU(3) breaking of 
the octet of (ir,K,ri) and can not be taken into account, a priori, in our way of fixing Dl(s) 



making use of lowest order xPT [16] and the exchange of resonances given by [24]. 

The fit will be done for the following experimental data: the elastic S-wave tttt phase 
shifts with 1=0, <5??, the KK -> vrvr, l=L=0 phase shifts, <Jg, the l=L=0 with 

r7oo the inelasticity in that channel, the elastic S-wave, 1=1/2 Kti phase shifts, 5h and 
a distribution of events around the mass of the ao(980) resonance, corresponding to the 



central production of ttttt] in 300 GeV pp collisions 143] for the 1=1, L=0 channel. 
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Because results coming from different experiments analyses are not compatible, we have 
taken as central value for each energy below y/s = 1 GeV the mean between the different 
experimental results jl^, 44]. For \/s > 1 GeV, the mean value comes from [44, [IjJ. In 
both cases the error is the maximum between the experimental one and the largest distance 
of the experimental values to the mean one. This procedure will be the one adopted, when 
needed, for the rest of the experimental magnitudes included in the fit. 
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For this quantity there are two sets of data below yfs = 1.2 GeV. Higher in energy both 
sets converge. One group will be represented by jlT]] and the other one by pq, 49]. The 
experimental results from 47] are larger than the data of the other works [48, ^] below 1.2 
GeV. We will distinguish between both cases when doing a fit referring it as high/below 
respectively. The change in the value of the fitted parameters will be very small when 
changing from one set of data to another, so that, this experimental ambiguity will not be 
relevant for our final values. We will average the experimental data of the second set of 
works for y/s < 1.2 GeV in the way explained above. When y/s > 1.2 GeV the average will 
be done between all the quoted analyses [47, |4^, 49]. 
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There are a series of analyses and experiments about the inelastic cross section -kit — ► 
KK, which agree between each other in the values for 1 . We have taken the data from 
1 47, [4^] as representative for such situation. 

The quantity 1 ?°° has been used instead of the inelasticity, r/oo, because the former is 



much better measured and all the experiments [47, 48, |4{|] agree on that quantity 



X2 " 



We distinguish between the more recent experiment |5C] and the older results [34, 
We have averaged the data from the last analyses up to y/s = 1 GeV. Above this energy, in 
the latter group of experimental works, only |35[] offers data. The statistical errors in this 
latter experiment are very small. We have enlarged them at the level of those in the most 
recent experiment |Hj which would make the different experiments compatible. Thus, the 
final points used in the fit for this magnitude will be the ones from |30| and the average 
between [34, [35], 51] as described above. 



1=1, L=0 data. 

The experimental data is very scarce for this channel. 



We will take a distribution of 

We 



events corresponding to the central production of mrr] in 300 GeV pp collisions [42] 
will study the data points around the mass of the ao(980), where one could think that the 
energy dependence will be dominated by the exchange of that resonance. We add in an 
incoherent way with respect the ao(980) resonance, the same background as in |f42]| . The 
ao(980) contribution is parametrized as: 



dN 

d-EUr, 



■ Np. 



7T7J I 



-CM 12 



(71) 
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with pnjj the three momentum of the ttt] state in the CM corresponding to a total energy 
E cm and N is just a normalization constant. 



The fit. 

Let us now discuss the fits obtained when including: (1) two scalar nonets with masses 
around 1 and 1.4 GeV or (2) only one nonet with a mass around 1.4 GeV. Of course, the 
final value for the tree level, 'bare', masses of the octet and singlet will be given by the fit. 
The fits have been done using the MINUIT minimization program. The output value for 
the parameters are written with the same precision as given by MINUIT. 

The fit that results when two nonets are included and also with the high 5™ data is: 



First Nonet (MeV) 

c d = 1.7997 

c m = 0.66534 

M 8 = 1003.41 ± 600 

c d = 20.988 

c m = 8.4867 

Mi = 1031.91 



Second Nonet (MeV) 



19.512 
= 19.612 
= 1379.31 

0.32994 
= -2.7183 
= 1000 ± 600 



a 



SL 



= -.72225 
9.2195 MeV" 



(72) 



Xd.o.f 



0.9688 
188 points 



A very striking aspect appears when observing the value of the parameters given in 
eq. (f72|). The value of the constants c<j, c m and c! d and c' m are, at least, one order of 
magnitude smaller than c d , c m and c' d , c' m respectively. This makes that the first octet 
and second singlet are phenomenologically irrelevant. Note that their masses are essentially 
undetermined. This is shown by the need to increase them by 600 MeV in order to make 
the Xdof t° increase by 0.5 units. In this way, they do not originate or participate in the 
poles corresponding to the physical resonances mentioned at the beginning of the section. 
They only give rise to poles very close to the real axis, with a width of only a few MeV. 
These poles manifest themselves as very narrow peaks in the partial waves, which are not 
observed by experiment. 

From the latter discussion we will just introduce a scalar nonet. The resulting values, 
after a new fit to the data, are: 



High 5°° 



Nonet (MeV) 
c d = 19.113±|f 
c m = 15.110 ±30 
M 8 = 1390.31 ± 20 
c d = 20.918±J;g 
c m = 1 0.567 -.^ 
Mi = 1021.11 



-10 



-20 



a 
N 



SL 



= -.75110 ±0.2 
9.3922 ± 4.5 MeV" 



Xd.o.f 



(73) 



1.066 
188 points 



Low 5^ 
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a 
Af 



SL 



= -.74246 
9.429 MeV" 



Xd.o.f 



(74) 



1.21 
196 points 



Nonet (MeV) 
c d = 19.183 
c m = 15.248 
M 8 = 1390 
c d = 20.941 
c m = 10.641 
Mi = 1021 

We have also shown the statistical errors for the parameters of the high 6®® fit obtained 
by increasing the Xdof by one un it, in order to appreciate the precision in the value of the 
parameters given by the last fits. The large error on c m is because this constant, as can 
be seen from eqs. (|66|), ( |68| ) and (|70|), enters through the multiplication of squared masses 
of the lightest pseudoscalars which are much smaller than s — M|, around the resonance 
region of the octet. Thus, its influence in the final value of the amplitudes is very small. 
This also happens to c m , although to a lower extension because Mi < Ms. One can see, 
comparing the last two fits, that the variation in the value of the parameters is very small 
when changing from one set of data to the other. The resulting fit for the high 5^ data 
is shown in Figs. 4-8. The results obtained before also favor the high solution for the 
phase shifts because its corresponding Xdo f * s sma U er than the one for the low solution. 

400 



300 




400 



800 1 000 

E cm [MeV] 



Figure 4: Elastic isoscalar tttt phase shifts, 5^. The circles correspond to the average of 
[43, 44, 45], as discussed in the 5®® subsection. We have also included the triangle points 



form [46 1 to have some data close to threshold, although these points have not been included 
in the fit because they are given without errors. 

This fit has 8 free parameters, 6 constants from the nonet^, a SL and the normalization 
constant M. In our former work [19], we were able to describe 5^1, d®® an d 1 up to 



If the singlet and the octet introduced form really a nonet is something we can not say. However, we 
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1000 1100 1200 1300 1400 
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Figure 5: S-wave KK — > tttt isoscalar phase shifts, <5^§ - The triangles points are from [47] 
circles correspond to the average of [48, 49 and squares to the one of 0, 48, 49 1. 



y^i < 1.2 GeV and a distribution of events around the ao(980) mass [52], using only 2 free 
parameters: a cut-off (which plays the role of the regularization scale \x and at the same 
time generates a concrete value for a SL ) and a corresponding normalization constant for the 
ao (980) event distribution. In fact, if we remove the resonance contributions in eqs. (BHI) 



and ( p7|) the formalism of ref. [19] follows from the present one. It might look surprising 



that a good fit to the data for -y/s < 1.2 GeV could be obtained in [19] with just one 



free parameter and a normalization constant for the mass distribution, while here one has 
needed 7 parameters, apart from the normalization constant. One reason is that now we 
have pushed the fit up to y/s = 1.4 GeV, while in |1| only data up yfs = 1.2 GeV was 
considered. The fact that new resonances appear around yfs =1.4 GeV has forced us to 
include an octet, which implies 3 new parameters, two couplings and a mass. However, the 
effect of this octet below y/s = 1.2 GeV is very small, hence only the singlet appearing with 
a mass around 1 GeV is relevant for energies below y/s =1.2 GeV. The present fit to the 
data has led us to the inclusion of this singlet resonance in T^° apart from the lowest order 
xPT Lagrangian, while in |TI| only the latter contribution was considered. The reason that 
forced us to include now the mentioned singlet is the consideration of the r]r] channel in 1=0, 
which was omitted in [pi| ], and is not negligible above 1 GeV, as can be seen in its strong 
coupling to the /o(980) resonance that we will obtain below. The ijrj channel affects mostly 
the mag nitude Should one have taken the available data for 7700 instead of those for 

1 1100 , which are measured with better precision from the tttt — > KK inelastic cross section, 



4 

the effect of the rjr] channel would be masked by the large errors in 7700 



will denote the global contribution of the introduced octet plus the singlet by using the word nonet as a 
shortcome. In this way, we also follow the nomenclature of p4j, inspired in the U(3) symmetry which holds 
for N c — » 00 . 
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It is quite interesting to recall that an 1=0 elementary state around 1 GeV has been 



predicted from QCD inspired models |53|, 54] and has also been advocated in phenomeno- 



logical analyses |l], Such state could be associated with the preexisting singlet state 
that we need. 

Resonances. 



Let us now concentrate on the resonance content of the fit presented in eq. (73). The 



octet around 1.4 GeV gives rise to eight resonances which appear with masses very close to 



the physical ones, /o(1500), ao(1450) and 1Tq(1430) |30[]. Thus, the correlation between tree 
level resonances, poles and physical resonances is clear around yfs ~ 1.4 GeV. However, 
this correlation is not so clear around 1 GeV. This issue will be the object of the following 
discussions 0. 

From Figs. 4 and 7, one can easily see two resonances with masses around 1 GeV, the 
well known /q(980) and ao(980) resonances. The first one could be related to the singlet 
bare state with M\ = 1020 MeV, but for the second we have not bare resonances to associate 
with, because the tree level resonance was included with a mass around 1.4 GeV and has 
evolved to the physical ao(1450). The situation is even more complex, because we also 
find in our amplitudes other poles corresponding to the /o(400 — 1200) = a and to the 
-KTq(900) = k. In Table 1 the pole positions of the resonances in the second sheet^j are given 
and also the modulus of the residues corresponding to the resonance R and channel i, C/ 2 , 
given by 

|CfCfl = Hm \( s -s R )T i3 \ (75) 

J s^s R 

where sr is the complex pole for the resonance R. 

While for the /o(980) one has a preexisting tree level resonance with a mass of 1020 
MeV, for the other resonances present in Table 1 the situation is rather different. In fact, 
if we remove the tree level nonet contribution from eqs. (65), ( ]67| ) and (|69|) the ao(980), 



a and k poles still appear as can be seen in Table 2. For the /o(980), in such a situation, 
one has not a pole but a very strong cusp effect in the opening of the KK threshold. In 
fact, by varying a little the value of a SL one can regenerate also a pole for the /o(980) from 
this strong cusp effect. In Table 2 we have not given an absolute value for the coupling 
of the /o(980) to the KK channel because one has not a pole for the given value of a SL . 
However, the ratios between the different amplitudes are stable around the cusp position. 
As a result, the physical /o(980) will have two contributions: one from the bare singlet state 
with Mi = 1020 MeV and the other one coming from meson-meson scattering, particularly 
KK scattering, generated by the lowest order %PT Lagrangian. 

In eqs. (]65|), ( |67|) and (^) when the resonant tree level contributions are removed, only 
the lowest order, 0(p 2 ), xPT contributions remain. Thus, except for the contribution to 
the /o(980) coming from the bare singlet at 1 GeV, the poles present in Table 2 originate 
from a 'pure potential' scattering, following the nomenclature given in |jT^ ]. In this way, 
the source of the dynamics is the lowest order xPT amplitudes. The constant a SL can 

4 We do not give a detailed study for the resonances with masses around 1.4 GeV because we have not 
included channels which become increasingly important for energies above ~ 1.3 GeV as 4-7T in 1=0 or Kr{ 
for 1=1/2. This makes that the widths we obtain from the pole position of the former resonances are 
systematically smaller than the experimental ones j^]. Thus, a more detailed study, which included all the 
relevant channels for energies above 1.3 GeV, should be done in order to obtain a better determination of 
the parameters for this octet around 1.4 GeV. 

5 I sheet: Im p\ >0, Im p2 >0, Im pz >0; II sheet: Im p\ <0, Im p2 >0, Im pz >0 
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Table 1: Pole position and residues for the full amplitude. 



1 — 

V S a = 


445 + i 221 MeV 


J U 


= 987 + i 14 MeV 


C°~ = 

S7T7T 


4.26 GeV 


c /o_ 


3.63 GeV 


S7T7r 


0.254 


S7T7T 

r u - 
s lf if 


0.51 


^VV — 


0.036 


a/0 
SJ7J7 

7fo_ 


1.11 


A> = 


1053.13 + i 24 MeV 




= 779 + i 330 MeV 


^KK 


5.48 GeV 


^Ktt 


4.99 GeV 


7^0- - 

K 


0.70 


(~ K 

<" K 


0.62 



Table 2: Pole position and residues when the bare resonant contributions are removed 



At 


= 434 + i 244 MeV 


yfif 


= cusp effect 


S7T7r 


4.21 GeV 


/■fo 
^KK 




S7T7T 


0.301 


rfo 

S7T7T 


0.38 




0.033 


7n 


1.04 


A> = 


1081.95 + i 13.3 MeV 


V*K = 


770 + i 341 MeV 


A^O _ 

^KK 


5.98 GeV 


^Ktt — 


4.87 GeV 


S7T?7 

Kit 


0.74 


_ 

t K 
SXtt 


0.61 
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Figure 6: 1 7 00 with rjoo the I=L=0 S-wave inelasticity. Circles |47j], triangles |48] 



be interpreted from the need to give a 'range' to this potential so that the loop integrals 
converge. These meson-meson states are shown in Fig. 9 in the chiral limit, setting all the 
masses of the pseudoscalars to zero and all the fp = f, where P denotes any pseudoscalar 
meson tt, K or r]. We see in this last figure a degenerate octet for 1=0,1 and 1/2 with a 
mass around 500 MeV and a singlet in 1=0 with 400 MeV of mass. In both cases these 
meson-meson resonances are very broad. 

The situation is very different to that of the former studied vector channels where all the 
physical resonances, the p and K* , originate from the preexisting tree level resonances. We 
already saw, at the end of the last section, when comparing the S and P-wave nir scattering, 
that the 0(p 2 ) xPT amplitude is 6 times larger for L=0 than for L=l around the resonance 
energy region. This implies that n-loops in L=l, with the 0(p 2 ) xPT amplitudes at the 



vertices, will be suppressed by a factor 



1 



6 n+l 



with respect to L=0. The suppression of loops 



is expected from Large N c QCD and this is in fact what happens for the vector channels, 
but for the scalar ones unitarity is unexpectly large, giving rise to these meson-meson 
resonances. 

As can be seen from eq. (|39|), these meson- meson poles, without tree level resonant con- 
tributions, originate from the cancellation between the inverse of the 0(p 2 ) xPT amplitude 
and the go function. As a consequence, the following relation between the masses of those 
resonances with / results 



M 2 (x f 2 /g 



(76) 



since go is 0(1) and / is 0{N C ) [25|, these masses will grow as N c . Thus for N c — > oo 



these resonances will go to infinity. This movement can be followed by suppressing the go 
function by a factor r from 1 (physical situation) to (N c = oo). It is then observed how 
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the resonances in Table 2, without the preexisting resonant contributions, disappear going 
to infinity. 




E om [MeV] 

Figure 7: Distribution of events around the ao(980) mass corresponding to the central 
production -Kirrj in 300 GeV collisions ||42j|. The abscissa represents the nrj invariant mass, 



E cm . The dashed line represents the background introduced in the same reference. 




1000 1200 

E om [MeV] 



Figure 8: S-wave 1=1/2 Kir elastic phase shifts, <^{ 20 . The triangles correspond to the 
average, as described in <5*{ 20 subsection, of |34|, 35, 51 1. Circles correspond to (5 
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Figure 9: Chiral limit. From left to right, Figs. 9a,b,c respectively. In Fig. 9a the poles 
of T found in the unphysical sheet are shown for 1=0. Analogously for 1=1/2 and 1 in Figs. 
9b, c respectively. 



5 Estimations of the unphysical cut contribution from 
and the exchange of resonances. 

In this last section we estimate the influence of the unphysical cuts for the elastic tttt and 
Kit S-waves with 1=0 and 1/2 respectively. The unphysical cuts will be approximated by 



means of %PT supplied with the exchange of resonances [24] with spin <1 in the t- and 
u-channels. The loops are calculated from %PT at 0(p ) and the exchange of resonances 
in the crossed channels accounts for a resummation of counterterms up to an infinite order, 
in the way explained in section 2 after eq. (p6|). The result can be taken directly from ref. 



where the tttt — * tttt and Kit — * Kit amplitudes are calculated up to one loop including 
explicit resonance fields [24]. 

In order to extract, from ref. p3| ], the contribution of the unphysical cuts, which we 
design by Ti, e f t , we have made use of eqs. (3.2), (3.10), (3.13), (3.14) for the tttt scattering 
and of eqs. (3.6), (3.16), (3.19) and (3.20) for the Kit one. We calculate the loop contribu- 



tions at the same regularization scale than in [23], that is, [i = M p = 770 MeV, the same we 



have taken in this work. In the following, when we refer to an equation in the form (m.n), 



it should be understood that this equation is the corresponding one from ref. [23|. The 



work of [23] contains loops and the exchange of resonances in the s, t and u channels. The 
exchange of resonances in the s-channel is also present in our work where the masses and 



couplings of the scalar resonances, eqs. (73), were fitted to data. Obviously the loops and 
the exchange of resonances in crossed channels, absent in our work, go to TLeft- On the 
other hand, we must also include in T^Left a polynomial contribution of 0(p 4 ) because the 
loop functions used in [23|, Jpg, M r P q and Jpq, where P, Q are tt, K or 77, and our loop 
functions, (— )go(s), differ in a constant. This polynomial contributions can be interpreted 
as subtraction terms from a dispersion relation of TLeft- Since the loops are calculated 
with ImT at 0(p 4 ) one needs three subtractions, which fix the order of the subtraction 
polynomial. Let us explain first the tttt scattering. 

From eq. (3.2) one has the expression of the elastic tttt amplitude with 1=0 in terms of 
the amplitude A(s,t,u), eq. (3.10). Making use of eqs. (3.2) and (3.13) the contribution of 
the loops in the s-channel is given by 

" — ^ J™( s ) + -^jiJkk{ s ) + 7^^w( s ) (77) 
It is straightforward to see that the imaginary part of eq. ([77]) is the one required by 
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Table 3: Influence of the unpysical cuts for the I, L=0 tttt and 1=1/2, L=0 Ktt partial 
waves. The three first columns refer to tttt and the last three to Ktt. 



V s 
MeV 


T 

L Left 

|To,u| 

% 


L Left 
r i lOO 

0,11 

% 


V s 
MeV 


T 

L Left 

|To,u| 

% 


T 

-L Left 
0,11 

% 


276 


3.7 


4.8 


634 


7.1 


8.7 


376. 


3.5 


5.1 


684 


3.7 


4.7 


476 


4.1 


5.7 


734 


0.3 


0.4 


576 


5.7 


6. 


784 


-2.5 


-3.3 


676 


8.1 


6.1 


834 


-5.7 


-7.2 


776 


11.2 


5.6 










unitarity up to 0(p 4 ) for the 1=0 S-wave tttt elastic partial wave with pions, kaons and 
etas as intermediate states. The squared amplitudes in front of the loop functions are the 
lowest order %PT amplitudes since loops are calculated at 0(p ). This is the same kind of 
result we would obtain for the loop contributions in the s-channel from the expansion of 
the generalization of eq. (|39|) to coupled channels up to the order considered in eq. ([77|), 
after dividing eq. ( |77| ) by a global factor 2 to match with our normalization in eq. 
a = 2. However, as we discussed above we use (—)go(s)u instead of Jfj(s) in eq. 
order to evaluate the loop contributions in the s-channel. Hence, we must include in ^Left 
the following expression 

1 2f f (J^s) + go(s)ii) + ^(r KK (s) + g (s) 22 ) + ^(r vv ( s ) + g (s) 33 ) (78) 

The former contribution, together with the exchange of resonances and loops in the 
crossed channels, after projecting over the S-wave using eq. (p]) with a = 2, define T^Left m 
our approach. 

For the 1=1/2, L=0 Ktt partial wave one has essentially the same situation than for tttt. 
From eqs.(3.6) and (3.19) one calculates the contribution of loops, which we project over 
the S-wave. The loops in the s-channel give the corresponding result to eq. ( |77| ) for the Ktt 
1=1/2 S-wave. In this case, instead of having the loop function Jp P (s), one can write it in 
terms of Jkit( s ) an d Jx-q{ s )- After taking into account the difference between Jpq(s) and 
our loop functions, one obtains the analog result to eq. ([78]) for Ktt. This contribution, 
together with the projection over the S-wave of loops and the exchange of resonances (eqs. 
(3.6), (3.20)) in crossed channels, give ^Left- 

We have not considered the tadpole contributions coming from pseudoscalar loops with- 
out flux of energy and the coupling of scalar resonances to the vacuum ( Fig. 2.b of [p3| ) 
because they are reabsorbed into the residues and positions of the CDD poles and the 
subtraction constant ao, eq.(|j), which we have phenomenologically fixed. 

Including explicit resonance fields as done in [^] increases the range of safe applicability 
of Chiral Symmetry from y/s ~400 MeV, accomplished in xPT, up to \fs ~ 700 — 800 MeV, 
as can be seen in [^3| when comparing their results with the experimental data. 

The results which we obtain for the contribution of T^eft i n the range of energies of 
[23] are shown in Table 3. In the second and fifth columns we show, respectively, the ratio 
between Tieft and the absolute value of our calculated I, L=0 tttt and 1=1/2, L=0 Ktt 
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partial wave amplitudes up to \/s « 800 MeV. In Table 3 we also compare ^Left with the 
tree level amplitudes Tg^. This ratio is also significative because the procedure which we 
have followed to arrive to a unitarized amplitude from Tq° 11 would not be much affected by 
the addition of T^Left which is a small correction with respect to T§° 11 . We see that these 
ratios are rather small. Therefore, this supports our point of view of treating the left hand 
cut as a perturbation in the range of energies we have considered. 

It is worth mentioning that this smallness of the unphysical cuts, as shown in Table 3, is 
a consequence of a cancellation between the contributions to T^Left from the loops and the 
exchange of resonances in crossed channels. In fact, the individual contributions in the irn 
case, for energies around yfs = 600 MeV, are of the order of 15-20% with respect to Tq° 1:l . 

In a recent work, ref. [p^] , the authors also combine the N/D method with chiral 
symmetry studying the (I,L)=(0,0), (2,0) and (1,1) n7r partial wave amplitudes. However, 
in this work only elastic unitarity is considered and the calculations are done in the chiral 
limit (m n = 0). On the other hand, the left hand cut is approximated only by the exchange 
of the p plus a scalar resonance without including loops in the crossed channels. These loops, 
as we have seen in this section, cancel to a large extend the crossed resonance contributions 
for the S-wave 1=0 rrir scattering. 



6 Conclusions 

Making use of the N/D method, we have developed the most general structure that an 
elastic partial wave amplitude has when the unphysical cuts are neglected. After matching 
this result with lowest order, 0(p 2 ), xPT |l(| and with the exchange of resonances with spin 
<1, in a way consistent with chiral symmetry as given in ref. [24], we extend the formalism 



to handle also coupled channels. Then, tttt and i^7r(I=l/2) P-wave amplitudes are described 
up to y/s = 1.2 GeV. It is shown that these amplitudes can be given rather accurately in 
terms of m n , mx, f and the masses of the p and K* resonances, when restrictions coming 
from Large N c QCD and unitarity are considered, in the lines of what was observed in Q. 

Next, the scalar sector is studied and good agreement with experiment up to \/s = 
1.4 GeV is found. An octet and a singlet are included with masses around 1.4 and 1 
GeV respectively. The former originates the observed /o(1500), oo(1450) and Kq(1430) 
resonances, the latter an important contribution to the physical pole of the /o(980). Other 
poles appearing in our amplitudes, the ao(980), a, k and an important contribution to the 
final /o(980), originate from meson-meson scattering with the lowest order xPT amplitudes 
plus the constant a SL as dynamics source. This situation is very different from the one 
observed in the vector channels where tree level structures dominate the scattering process 
and a strong suppression of unitarity loops occurs, as indicated at the end of section 4. 
As a consequence, the present study supports that a concept like scalar meson dominance, 
analogous to the well known vector meson one, is not suited at the phenomenological level. 

In the last section we have made some estimations in order to investigate the influence 
of the unphysical cuts. The results obtained support our picture of treating the unphysical 
cuts in a perturbative way and then establishing the stability of our conclusions in sections 
3 and 4 against the corrections coming from cross symmetry. 



Acknowledgments 

We would like to acknowledge fruitful and basic discussions for the present work with 



31 



A. Pich. Discussions with S. Peris, J.R. Pelaez and A. Kataev are also acknowledged. This 
work was partially supported by DGICYT under contacts PB96-0753 and by the EEC-TMR 
Program- Contact No. ERBFMRX-CT98-0169. J. A. O. acknowledges financial support 
from the Generalitat Valenciana. 



A N/D in coupled channels 

In this appendix we make use of a matrix formalism to deal with several coupled channels. 
In analogy with the elastic case, eq. ([?]), let us define the matrix T' L as 

T' L (s)=p- L T L ( S )p- L (79) 

with p a diagonal matrix which elements are pij = piSij where pi is the modulus of the CM 

to? • , mk- ) , 
momentum of the channel i, Pi = 7= , with mu and rri2i the masses of the two 

mesons in channel i. 

From the beginning we neglect the unphysical cuts. As a consequence Ti(s)y will be 
proportional to pfpj- This makes that Ti(s)ij, apart from the right hand cut coming from 
unitarity (above the thresholds for channels i and j, s\ h and si, respectively), will have 
another cut for odd L between s\ h and s 3 th due to the square roots present in pi and pj. 
In this way, T' L will be free of this cut and will have only the right hand cut coming from 
unitarity. Thus it will satisfy 

Im T'jT 1 ^) = -p L p(s)p L = -p(s)p 2L (80) 
where p(s) is a diagonal matrix defined by 

p(s) = -^0(s) (81) 

with 9(s) another diagonal matrix such that 8(s)n=l above the threshold of channel i and 
below it. 

We write T' L as a quotient of two matrices, Nj, and making use of the coupled 



channel version of the N/D method |29] 



T' L = D'^N'l (82) 

We can always take free of poles and also containing all the zeros of T' L . In such a 
case N^ will be just a matrix of polynomials, we then write 

N' L = Qn-L-t (83) 

with Q n _£_i a matrix of polynomials of maximum degree n — L — 1. 
In this way, from eq. (BO) and (|82|) one has 



Im D' L ( S ) = -N' L (s)p(s)p 2L (84) 
and making a dispersion relation for T>' L one has 



32 



with P n _i a matrix of polynomials of maximum degree n — 1. 

Because is just a matrix of polynomials, it can be reabsorbed in D' L to give rise to 



a new T> L which will fulfill eq. (84) but with = 1. In this way 



L = ° L 
Nr. = 1 



-/ 



BL = * Jo (s'-s)(s'-s Q )^ +R ^ (86) 

with R(s) a matrix of rational functions whose poles will contain the zeros of T' L . This fact 

is in clear analogy with the role played by the CDD poles included in Section 2 for the 
elastic case. 
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